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Abstract 

Random Hermitian matrices are used to model complex systems without time-reversal invariance. 
Adding an external source to the model can have the effect of shifting some of the matrix eigenvalues, 
which corresponds to shifting some of the energy levels of the physical system. We consider the case 
when the n x n external source matrix has two distinct real eigenvalues: a with multiplicity r and zero 
with multiplicity n — r. For a Gaussian potential, it was shown by Peche [32] that when r is fixed or 
grows sufficiently slowly with n (a small-rank source), r eigenvalues are expected to exit the main bulk 
for a large enough. Furthermore, at the critical value of a when the outliers are at the edge of a band, 
the eigenvalues at the edge are described by the r-Airy kernel. We establish the universality of the r-Airy 
kernel for a general class of analytic potentials for r = C(n 7 ) for < 7 < 1/12. 
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1 Introduction 

Fix a Hermitian matrix A. Equip the space ofnxji Hcrmitian matrices M with the probability measure 

^dM) = _l e -"Tr(V(M)-AM) dM . ^ . = f e -nlV(V(M)- AM)^ 

z n J 



where dM is the entry- wise Lebesgue measure and the integration is over all Hermitian matrices. The 
eigenvalues of M represent the energy levels of a system without time-reversal invariance [35]. When the 
external field A is nonzero and V(M) = M 2 /2, this measure arises in the study of Hamiltonians that can 
be written as the sum of a random matrix and a deterministic source matrix [18]. 

When A = (no external source) and V(M) = M 2 /2, (1-1) describes the Gaussian Unitary Ensemble, 
or GUE. For reasonable choices of V(z) the spectrum tends to accumulate on fixed bands on the real axis. 
Introducing an external field A can have the effect of perturbing the expected position of the spectrum. For 
example, Aptekarev, Bleher, and Kuijlaars [16, 3, 17] studied the Gaussian case when the matrix A has two 
eigenvalues ±o, each of multiplicity n/2. When a is sufficiently small (the subcritical case), the eigenvalues 
of M accumulate with probability one on a single interval, just as when A = 0. As a increases the interval 
splits into two (the supercritical case) . There is a transitional value of a between these two cases (the critical 
case) where the local eigenvalue density near where the bands are about to split is described by the Pearcey 
process. See [31] and [14] for further studies of large-rank external field models, and [2] for recent results on 
the universality of the Pearcey process. 

We are interested instead in small-rank sources of the form 

A = diag(a, . . . , a, 0, . . . , 0) (1-2) 

r n—r 

assuming that r(n) = C(n 7 ), < 7 < 1/12. The ratio of r to n, which is asymptotically small, will be 
denoted as 

k '.= — = Oin^ 1 ). (1-3) 

The limiting distribution of the largest eigenvalue for such a small-rank external source in the Gaussian case, 
V(M) = M 2 /2, was studied by Pcche [32]. Again three distinct behaviors were observed. For a sufficiently 
close to zero, i.e. the subcritical case, the largest eigenvalue is expected to lie at the right band endpoint 
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and behave as the largest eigenvalue of an n x n GUE matrix. For a large enough, i.e. the supercritical 
case, r eigenvalues are expected to exit the bulk and be distributed as the eigenvalues of an r x r GUE 
matrix. In the transitional critical case, when the outliers lie very near the band endpoint, the distribution 
for the largest eigenvalue is an extension of the standard GUE Tracy- Widom function [34] when r = (see 
also [4, 5, 1]). These functions, denoted by F r (x), were first discovered by Baik, Ben-Arous, and Peche [4] 
in the context of sample covariance, or Wishart, matrices and were shown to be probability distributions 
by Baik [5]. Adlcr, Delcpinc, and van Moerbeke [1] showed that these distributions also appear when n 
non-intersecting Brownian motions start from x — when t = with n — r conditioned to end at x = 
when t = 1 and r conditioned to end at x — a when t = 1. The walkers all start out in a single group, but at 
a critical time depending on a, a group of r walkers separates from the main bulk. At this critical time the 
walkers on the edge where separation is about to occur follow the r-Airy process, which is connected with 
F r (x). See [5] for other processes in which the functions F r (x) arise. 

In this paper we extend Peche's result in the critical case to more general functions U(M). Our specific 
assumptions are listed in Section 1.3, but we essentially assume U(M) is a generic analytic potential with 
sufficient growth at infinity. This establishes a new universality class of matrix ensembles with the local 
eigenvalue density near the critical point described by the r-Airy process. The universality of the supercritical 
and subcritical cases has been considered separately [9]. 

In the case of rank one perturbation (i.e. r = 1), the recent paper [7] by Baik and Wang has described 
the limiting distribution of the largest eigenvalue for all the possible cases including the critical case that we 
consider here. 

1.1 The kernel and its connection to multiple orthogonal polynomials 

Let p m (Xi, . . . , A TO ) be the probability density that the nx n matrix M chosen using (Tl) has eigen- 
values {Ai,...,A m } (here m < n). Then, when the Aj are distinct, the m-point correlation function is 
( w "| ra )i Pm(Ai, • • ■ , A m ). Brezin and Hikami [18, 19, 20, 21] showed that in the Gaussian case, the m-point 
correlation functions can all be expressed in terms of a single kernel K(x, y): 

Til 

— T p m (Ai, . . . , A m ) = det(K(Xi, \j))i,j=i,..., m . (1-4) 

(n — my. 

Zinn- Justin [36, 37] extended this result to the case of more general V(M). We will find the leading term in 
the large-n asymptotic expansion of the kernel in the critical regime near the critical endpoint. 

Bleher and Kuijlaars [15] showed that the kernel can be written in terms of multiple orthogonal poly- 
nomials. Furthermore, these multiple orthogonal polynomials can be written in terms of the solution to a 
certain Riemann-Hilbert problem. Specifically, suppose Y(z) is a 3 x 3 matrix-valued function of the complex 
variable z satisfying 

'Y(z) is analytic for z ^ R, 

(I e -nV(x) e -n(V(x)-ax)\ 
1 forxeM, 

1 ) (1-5) 

(z n ( 

Z ( " r ' ! /r- • - X . 
y Z 



Y(z) = (l + G(i)) 
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Here Y±(x) := lim^o Y(x± is) denote the non-tangential limits of Y(z) as z approaches the real axis from 
the upper and lower half-planes. Whenever posing a Riemann-Hilbert problem we assume (unless otherwise 
stated) that the solution has uniformly Holder continuous boundary values with any exponent p € (0, 1] 
along the jump contour when approached from cither side. Under our assumption (iv) in Section 1.3, the 
unique solution Y(z) can be written explicitly in terms of multiple orthogonal polynomials of the second 
kind (see [16], Section 2). In the case of two distinct eigenvalues a and 0, which is our case, the kernel 
K n (x,y) may be written in terms of the function Y(z) as 

e -in(V(s)+V(i,)) A\ 

K n (x,y)= 2m{x _ y) (0 1 eUay ) Y(y) Y(x) I j . (1-6) 

To analyze the asymptotic behavior of Y we will use the standard nonlinear steepest descent method for 
Riemann-Hilbert problems, as well as certain ideas introduced by Bertola and Lee [10, 11] to study the first 
finitely many eigenvalues in the birth of a new spectral band for the random Hcrmitian matrix model without 
source. 

A potential alternate method for establishing universality for finite r would be to use Baik's result [6] 
writing the kernel K n (x, y) in terms of the standard (not multiple) orthogonal polynomials. The rank of the 
matrices in this alternate expression grows with r, whereas the size of the Riemann-Hilbert problem (1-5) 
grows with the number of distinct eigenvalues. As such, for growing r it is more convenient to analyze the 
Riemann-Hilbert problem for multiple orthogonal polynomials. 

1.2 Definition of the critical regime 

We recall the setting of our work [9]. Let g(z) be the (/-function associated with the orthogonal polynomials 
with potential V(z) (see, for instance, [23] or [25]). It may be written as 

q(z; k) := - — — - / log(z - s)p min (s; n)ds = log(z) + O ( - ) , (1-7) 

where Pminds = p m - m (s; n)ds is the unique measure minimizing the functional 

T[p]:= I V(s)p(s)ds- [ [ p( S )p(s')log\s-s'\dsd S ' , [p(s)ds = l-^. (1-8) 

JR JR JR JR * 

Here n is a small parameter which is identified with the ratio k = - . The variational equations are equivalent 
to the statement that [33] there exists a real constant l\ = li(n) such that 

Re[V(x)-(2- K )Q(x;K)-t 1 (K)]\ ^ ° on M \ suppp min , 
L ' x ;i y = onsupp/9 min . 

We shall denote 

g{z) :=Q(z;0) , ii := ^i(O). (1-10) 

Our first assumption will be 

Assumption 1.1. The unperturbed (k = 0) variational problem is regular in the sense of [30] which means 
that the inequalities in (1-9) are strict and the behavior of V(x) — 2g — li at any boundary point x = £ of the 
support of p is asymptotic to oc (x — £)? (approaching £ from the complement of the support). 
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It has been shown in [30] at Theorem 1.3 1 that, for real-analytic V(x), 

i) If V(x) is regular for k = then V{x) is still regular for n small enough, and 

ii) The locations of the spectral edges (the a's and /3's) are real-analytic functions of k such that the bands 

of the support of the equilibrium measure stay separated as k ranges in a small open set around k = 0. 

In addition it is also known that the support (under the real-analyticity assumption) consists of a finite 
union of bounded intervals; we will denote the support of the density p min by (see Figure 2) 

suppp min = ^(J U [oi(k),I3(k)], ^ ^ 

«i(k) < /?i(k) < a 2 {n) < (3 2 (k) < • • • < a g (n) < /3 9 (k) < a(n) < /3(k). 

We will consider the unperturbed density to be the solution of the above variational problem with n = 0; in 
this case the reference to k will be tacitly suppressed, and so otj — ctj(K,)\ Q , etc. Define for the unperturbed 
(k = 0) problem the following quantities: 



P2(Z) 

Ps(z) 



= -V(z) + 2g(z) + h, (1-12) 
= -V{z)+az + g{z) + h-h, (1-13) 
= -P 1 (z) + P 2 (z) =az-g{z)-h. (1-14) 



Note that V(P) = \g{$) - h and hence P 2 (/3) = P 3 (/3); we choose l 3 such that P 2 (/3) = P 3 {/3) = 0. 

It is also known that Reg(z) is a continuous function on R and harmonic (and convex) on the complement 
of the support (up to a sign it is also known as the logarithmic potential in potential theory). 

Definition 1.1. Define a c to be the (unique) value of a so that P 2 (f3) = (here n = 0). 

The uniqueness is promptly seen because P 2 ((3) — —V'(f3) + a + g'((3); in fact the effective potential Pi 
is known [24] to satisfy 

P[(z) = 0(z-p)K (1-15) 
In particular, P{(/3) — and hence P 2 ((3) = a — g'((3) = a — ^V'((3). Thus the critical value of a is given by 

a c = g'(P)= 1 2 V'(f3). (1-16) 

We also recall that for regular potentials the behavior of (a suitable branch of) the function P\{z) near 
any of the endpoints of the interval of support is 

P 1 (z) = -C(z-p)i(l + 0(z-p)) (1-17) 

for some constant C. For the point ft — sup [supp p m i n ] one can also prove that C > 0; this allows us to 
introduce the scaling coordinate £ near z — /3 via the definition 

2 

( = ((z;n):={-—P 1 (z)) =n*a(z- 0){1 + O(z - 0)) , a > 0. (1-18) 



-'^In the theorem, the changing parameter is essentially k after rcscaling. 
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We now define the critical and near-crtical regimes. A more extensive context for these definition can be 
found in [9]. For completeness we also define the supercritical, subcritical, and jumping outlier regimes. The 
supercritical and subcritical regimes are dealt with separately in [9]; we plan to consider the (non-generic) 
jumping outlier regime in a future work. 

Definition 1.2. The matrix model specified by (1-1) is in the critical regime if a = a c and P2{x) < ^(/S) 
for x > j3. The scaling regime of a — a c = O (n -1 / 3 ) will be called near-critical. We define the exploration 
parameter r by 

t := lim n 1/3 (a- a c )/d, (1-19) 

n— >oo 

where c\ is the positive constant defined at (1-18). 

We define, for < a < a c , b* to be the unique point on the real axis greater than /3 such that P^((3) = 0. 
For a > a c we can choose b* := (5. 

Definition 1.3. The model is in the supercritical regime if P 2 has a unique global maximum on {x > 
max{/3, b*}} at a point x = a* e R and any of the three conditions below is satisfied: 

• a > a c . 

• a = a c and Pz{(3) = Pi{fi) < Pi{x) for some x > (3. 

• < a < a c and Ps(b*) < P2(x) for some x > b* . 

Note that a* is always greater than (3 and b* . If the global maximum on [max{/3, b*}, 00) is attained at several 
distinct points then we will say that we are in the jumping outlier regime. 

Definition 1.4. The matrix model specified by (1-1) is in the subcritical regime if a < a c and P2(x) < 
P 3 {b*) for allx> b*. 

1.3 Assumptions and results 

We will make the following assumptions on a, A, and V(z): 

(i) a > 0. 

(ii) A is a small-rank external source of the form (1-2) with r = C(n 7 ), < 7 < 1/12. 

(iii) V(z) is real analytic and regular in the sense of [30]. 

(iv) hm - — — — — = 00 and lim - — — — — = 00. 
|z|-kx> log(l + z 2 ) \z\^oo log(l + z 2 ) 

Regarding assumption (i), the case when a < is equivalent by sending a — > —a and V(z) — > V(—z). 
As for assumption (ii), in the general case when A has m > 2 distinct eigenvalues the kernel can be 
written in terms of multiple orthogonal polynomials associated to an (m + 1) x (m + 1) Riemann-Hilbert 
problem, which is beyond the scope of this paper. The assumption of analyticity in (iii) allows us to use 
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the nonlinear steepest-descent method for Riemann-Hilbert problems. The assumption of regularity ensures 
that the equilibrium measure of V(z) has square-root decay at each band endpoint (so that we can use Airy 
parametrices) and that these endpoints are analytic functions of k near k = 0. 

Next, (iv) guarantees the existence of the multiple orthogonal polynomials needed to ensure the Riemann- 
Hilbert problem has a solution. We note that the allowed V(z) include any convex V(z) (see the introduction 
of [9]). 

We compute the large- n behavior of the kernel function (1-6) in the critical regime. We explicitly compute 
the kernel in a neighborhood of /?. In the remaining portions of the complex plane, our result is that the 
kernel function converges to the kernel for the classical orthogonal polynomial problem with respect to V(x). 
That is, away from /3 the standard universality classes apply (i.e. the sine kernel in the bulk of the spectrum 
and Airy kernels at the other edges). Our main result is: 



Theorem 1. Suppose V(z) and a satisfy conditions (i)-(iv). Let Cx,Cy oe fixed in some bounded set. Also 
let c\ be the constant appearing in (1-18). Then for large n, and for a positive integer r such that r = 0(n J ) 
with < 7 < 1/12, 



a. Cy+S 
2 ; P ' 1 



C\n 



2/3 



(2m) 



;\2 ~ 



ds / dt 



(t + T ) r e s(" 3 - f3 >+C- f -C» B / / 1 

{s + T) r t-s \ + \nV 3 



(1-20) 



where the oriented contours C and C are given in Figure 1 and S is a quantity independent of £ x , Cy and of 
the form 



S := Cx/Skti 3 = £>(n 7 a) 



0:-- 



dj8(«) 
dn 



(1-21) 



t=0 



Remark 1. By dropping the drift term S in (1-20) we would simply deteriorate the error estimate to OirC 1 3 ) 
which is however still vanishing since 7 < ^ . 

The constant /3 admits an explicit integral representation for an arbitrary real-analytic potential V(z) 
but it is a bit complicated when the equilibrium measure is supported on multiple intervals. In the simplest 
case where the support of the equilibrium measure consists of a single interval [a, 0\ then we have 

1 , . , (1-22) 

la \ X V'(z)dz ' v ' 

where the contour of integration is a simple closed contour surrounding the support [a, /3] in the complex 
plane. We will not be using in any way the explicit form of /?, except the fact that it is a well-defined 
quantity due to the smoothness of (3(k) guaranteed by the already cited Kuijlaars' Theorem 1.3 in [30]. 

We show in Section 6.3 that our kernel is the same as the one found for nonintersecting Brownian walkers 
in [1]. 

Acknowledgments. The authors would like to thank Jinho Baik, Ken McLaughlin, and Dong Wang for 
several illuminating discussions. M. Bertola was supported by NSERC. R. Buckingham was supported by 
the Charles Phelps Taft Research Foundation. V. Pierce was supported by NSF grant DMS-0806219. 
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Figure 1: The contours C and C used in the r-Airy Kernel. The straight lines cross at — t. 

2 The perturbed equilibrium measure and the local coordinate £ 



To describe a growing number of outliers, we define a perturbed equilibrium measure problem with parameter 
of perturbation k = r/n = C(n 7_1 ) for < 7 < 1/12. Using the g-function (1-7) we define 



Vi (x; k) 


:= -V(x) + 


(2 - k)q(x;k) +£i, 






V 2 {x; k) 


:= -V(x) + 


ax + -(2 — k)q(x; k) + 


^l og (x-^) + 4-4, 


(2-1) 


V s (x; k) 


:= -V1+V2 


= ax--(2 - k)q(x; k)- 


h ^log(x-^)-4. 





The other constant £3 will be defined by (2-7) after we define the locally holomorphic function h{z). 
We also note that we have 

\imV' 3 (x) = P;(x) 1 j = 1,2,3, (2-2) 

uniformly over compact subsets of K \ {/3}. The Pj(a;)'s (j — 1,2,3) are defined by the same equations 
as the 7-ys except that n is set to zero (hence all the log terms are dropped) and that g(x) is replaced by 
g(x) = lim K _>.o q{x). Such convergence is uniform (outside a finite disk around /? and inside a compact set) 
according to [30]. 

Since V is regular for k = 0, it will still remain regular for small values of k. Hence there exists a 
holomorphic function £(z) in a finite disk around (3 such that 

-nV 1 (z: K ) = ^((z; K ) 3 / 2 , k)=0, (2-3) 
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and thus 

C = n 2 ^ Cl ( K )(z - /3(«))(1 + 0(z - /?(«))) (2-4) 

for some constant Ci(k) > 0. 

Previously, we have defined a c at k = 0. To describe the effect of growing r we introduce a more exact 
definition a c {n). Note that 

h(z; k) := \v x {z- k)+V 3 (z; K )-^log( = - l -V{z) + az + ^ log + ^ - £ 3 (n) (2-5) 

is a locally holomorphic function at z = because V(z) is real analytic. 

Definition 2.1 {a c {n) and £s(k)). For n > 0, a c (n) is such that h'{P(n); k) = at a = a c (n), i.e. 
FFe a/so define i 3 {n) smc/i t/iai 

W«);«) = i(-Iogr(r)-log(^)) I (2-7) 
where T(r) is Euler's Gamma function and 

Vo(r) := J™ Ai(t) (e Tt '" + e wrt + e Tt ) dt , « := cxp . (2-8) 

One /mds explicitly 

««) := i (>o g r W +1 „ g - + - ^ + ^ lo g (^) ^ (2-9) 

We observe that if r grows with n then Re(/i(/?(/t); n)) < for sufficiently large n, which will be used in 
the proof of Proposition 5.5. One can verify that a c (n) converges to a c := a c (0) (see Definition (1.1)) when 
k — ^ at the rate 

a c ( K ) = a c + <D(n) =a c + 0(n*- 1 ) . (2-10) 

Since for k — we have a c — \V'{(3) — </(/?) > 0, for sufficiently small k we still have a c (n) > 0. 

From Definition 2.1, we have h(z) — h((3(n); k)+0((z—(3(k)) 2 ) at a = a c (n) (because h'{(3(n); k) — 0). For 
other values of a, since only the term az in h(z; n) depends on a, we get h(z; n) — 0((z— /3(/-c)) 2 ) + (a— a c {n))z. 

For the subsequent exposition, we redefine r in a way that is compatible with the earlier Definition 1.2. 

Definition 2.2 (replacing Definition 1.2). 

t ■.= n 1/3 (a-a c {K))/c 1 (n). (2-11) 
From (2-4) and the above definition, we get 

h(z; k) = h((3( K ); k) + -C + ((z - f3( K )) 2 ) = h(f3( K ); k) + ^r(z - /?(«)) + 0((z- f3(n)) 2 ) . (2-12) 
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3 Initial analysis of the Riemann-Hilbert problem: the global 
parametrix 



We define the contour L as the positively-oriented circle centered at a± (the leftmost edge of the spectrum) 
and passing through (3 (the rightmost edge of the spectrum). We choose the circle L large enough so that 
V2 is negative on the real axis to the left of L. Until further notice the dependence of the various quantities 
on k (i.e. a c , f3, g, etc.) will be understood throughout. We have: 

Lemma 3.1. For sufficiently small k and a — a c , the function Re(— V^{z) + 4p log(z — /?)) increases as one 
follows L in either direction starting from f3 (i.e. through the upper half-plane or the lower half-plane) . 

Proof. From the definition (2-1), we have —Vz{z) + 4^ log(z — )3) = —az + |(2 — k)q(z) — I. It is obvious 
that Re(— az) increases when a = a c > 0. It is also simple to see that Re (|(2 — k)q(z)) increases along the 
referred contour because Relog(z — x) increases along the contour for any x in [ai,/?]. □ 

We now open up lenses around each of the bands in the standard way (as in the analysis of the orthogonal 
polynomials associated to V(z)). We introduce the following open regions (see Figure 2): 



ttj , j = 1, ...,g: The area in the upper half-plane (+) or lower half-plane (— ) between the band [c(j,(3j] 



± 

and its appropriate adjacent lens. 

^main- The area m the upper half-plane (+) or lower half-plane (— ) between the band [a,0\ and the 
appropriate adjacent lens. 

ftf: The area in the upper half-plane (+) or lower half-plane (— ) inside the contour L but outside the 
lenses. 

fi^ ut : The part of the upper half-plane (+) or lower half-plane (— ) outside the contour L. 
We also define 




aens := n+ ns u n lcns , n L -.= n+un L , out :=n+ t un o - ut . (3-1) 

With these definitions of regions and contours, we define W(z) in each region by 

{ e-i Vl 

J(z) I et Pl I e -f«l°s(*-« 5 (3_2) 



' e ir*i 

e-* £l ) (3-3) 

et^-'O, 




-fV 

























and 




JC):= { I, 

z € fiout z E FIl z E 




(3-4) 



lens 
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Then W(z) satisfies the following jump conditions (see Figure 2 for the orientation of the contours): 

z £ L = dfl out n <9f2 L , 



W+{z) = W_(z) x < 








o \ 


n 
u 


1 -e^W , 








1 / 




1 


0\ 




nT'i(z) 


1 , 







1 







(-l) r 


(- 


-l) r+1 










1 


e" 


incr(z) ^_ 


ly (-i 







(- 













e nV 2 {z) 





1 














z e <9f2i ns n dfi L , 

z& (Uf=iK-.ft])u[a,/3], 



zeSfi+nSJl,;, 



zeafi+ ut n9n- t) 



(3-5) 



and the boundary condition 



W(z) = 1 + 



oo. 



(3-6) 



Here we have defined ct(x) by cr(x) := f*^ p m in{s; n)ds which is some constant on each gap. 

We will show below that the above Riemann-Hilbert problem is exponentially close to a simpler one away 
from the turning points such as (3. To be more precise, let us define a shrinking disk centered at (3 by 



< n 



-1/2 j 



Here we have chosen the diameter of the disk D c so that (refer to (2-12)) 

n(h{z)-h{P)- T -C) =nO((z-pf) 
is uniformly bounded on D c as n — > oo. From the expansion (2-4) and (3-7) we note 

C = ©(n 1 / 6 ) on dD c . 



(3-7) 



(3-8) 



(3-9) 



To complete the error analysis we will also need a fixed-size disk around /3. Fix a small <5 > and define 
D(/3, (5) to be the finite disk centered at /3 with fixed radius <5. For n large enough, we have D c C D(/3, S). 
We will show in Section 5 that, for z outside D c and a finite distance away from the other turning points 
(the Qj's, /3j's, and a), the Riemann-Hilbert problem for W(z) is exponentially close to that for the outer- 
parametrix ^ as n — >• oo with the proper choices of l\ and £3. The necessary data on the effective potentials 
V\ (z), T'2{z), and V${z) is contained in Lemma 5.5. Thus we propose the following model Riemann-Hilbert 
problem for the outer parametrix. 
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Definition 3.2. In the critical case, the outer parametrix *&(z) is defined as the solution of the Riemann- 
Hilbert problem 

*+(*) = *_(*) I (-1T+ 1 1 , z e (Uf=iK-,ft]) U [a,/?], 

/(-l) r e- i,l<T ( z ) 0\ 

*+(«) = *-(*)( (-l)'V-(z) Oj , «eR\((Uj=i[aj ) ft])u[a ) ^]), (3-10) 

*(z) = O ((z - *T 1/4 ) , z^*, * = a j ,P j ,a,p. 

*(«)=I + 0(i), Z^OC. 

TTie function a(z) is defined below the equation (3-6). 

This Ricmann-Hilbcrt problem is essentially 2x2. The unique solution is given in [25] Lemma 4.3, or 
in [13] Lemma 4.6 in a slightly more generalized setup. We refer the interested reader to those papers since 
this is not essential here. In the subsequent analysis we will only need the following information. 



Lemma 3.3. Define 

(( c \" 1/4 

S(C;n) := 



( r»V3 ) 

(^) 

V o 01/ 



1/4 

[-^r] 



'(-l) r 

U ( 1 I ; U := ( -\ f ) . (3-11) 
(-l) r 




Then for z € D C7 f/iere is a unique holomorphic 3x3 matrix-valued function H( )(z;n) determinant 
one such that 

M>(z) = H (0) (z)S(C) (3-12) 
as £ — > oo. In addition, H( )(z;n) /ias a Izmit asn^oo and 

H (0) (z;n)- limH (0) («;n)=O(K). (3-13) 

Proof. A direct check shows that S(£;n) has the same jumps as \&(z) in a neighborhood of z = /3. Thus 
the ratio VPS 1 has no jump discontinuities inside D c and hence may have at most an isolated singularity at 
z = (3. Furthermore, this product has at worst a square-root singularity at /3 (coming from the product of 
the quarter-root singularities in * and S). In the absence of a branch cut, this means H( ) is holomorphic. 
Finally, (3-13) follows from the definition of Q(z) in (2-3) and the dependence of V\ on k. □ 

4 The local parametrix near z = (3 

We begin this section by expressing the Riemann-Hilbert problem satisfied by W inside D c in terms of the 
local coordinate Q. Zooming in on D c , the contours are shown in Figure 3. There we collapse the global 
contours L and a part of <9f2 ma in into T 2 and IV The regions II and III are parts of the region fi ma i n - 
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Figure 3: Contours where A r has jumps. The location of the center can be in any finite domain. 
Figure 4: Contours to define the generalized Airy functions. The center is located at — r. 



Using the identities below, 



v 2 {z) = \v 1 {z) + h{z) + y 1o sC, 

Vs{z) = -\vi{z)+h{z)+^\ogt, 



we see that W(z) satisfies the following Riemann-Hilbert problem inside the disk D c : 



'1 e~t^ /2 ^ e -§C 3/2 +^W N 
1 
,0 1 



1 







_ e fC^ i _ C f e §C 3/2 +nM*) 



4^3/2 





W+(z) = W_(z) < 



(-l) r ON 
(-l) r+1 
1, 



1 

4^3/2 

es^ 1 — Q 2 e 









f p §C 3 / 2 +r>fc(*) 



zeT 2 , 



z e r 3 , 



z g r 4 . 



We construct below a local parametrix that solves a Riemann-Hilbert problem similar to the above. 
4.1 r-Airy parametrix 

The rth derivative of the standard Airy function admits the contour integral representation 

|fAi(C) = ^- / fe^'Ht, r = 0,1,2,..., 



14 



where the contour C\ is shown in Figure 4. Extending the standard Airy function, we will need the following 
generalized Airy functions. 

Definition 4.1. Let us define the following generalized Airy functions corresponding to each contour in 
Figure 4- 

Aig(C;r):=^ (t + rYe^^dt, (4-5) 
where m — 1, 6 indicates the contour depicted at Figure 4, and r is a non-negative integer. 

Using the generalized Airy functions, we shall construct a matrix A r (C) satisfying the jump condition 
(see Figure 3) below. 

Cer 1 (er 2 (er 3 (er 4 

(4-6) 






Ar(C)+=Ar(Q-X 

Theorem 2. For a positive integer r, the following definition satisfies the jump condition in (4-6). 



MO 



4 r) (0 4 r) (0 



d ( v[ r \o d c vi r \Q 9 c 4 r) (0 
U r_1) (0 vt 1} (o vt 1] (o 










)(: 


1 


i) 










where 



'Ai^(C), (el and II, 



Ai^ ; (C), (£ III and IV, 



-1 for II, 
1 for III, (4-7) 
for I and IV, 

a4?(0. Ce/ 

Ai«(C), (ell and III, 

[Ai«(0, Ceiv. 

(4-8) 



Proof. By applying the Sokhotskyi-Plemelj formula, one can satisfy the jump condition (4-6) by defining 

v[ r \() :=Aig(C;r), 

v[ r \s)e TS 



4 r) (C) : 



er T < f v { ;\t)e Tt dt 



(4-9) 



t-C 



2iT7 



v { 2 r \t)e Tt dt 



Here we note that Ai^(C; r) is holomorphic (therefore having no jump). The next step is to verify that the 
definitions in (4-9) are equivalent to the advocated form in (4-7) and (4-8). 



First let us consider v^\(). 




ds 1 

s — £ 2in 



dt + 

t — ioo 



C f~ T f rioo — r 

ds dt+ ds dt 

JSC— J — T—ioo J SC+ J—t 



dt)(t + rYe s ^- t3 ^-A^(C;r) 
(t + T) r eA t+ ^- t3 / 3 



(4-10) 



(t + rYe^/'dt^A^^r), (eM+, 
l -J c (t + rye^- t3 / 3 dt = A^(C;r), ( e H~, 



+ — f \t + T) r e <:t ' t ^ 3 dt (4-11) 

2«7T Joo 



(4-12) 
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where, at the second equality, sc- and sc+ are respectively the completions of the contour R by adding the 
infinite half-circle through the lower (— ) half-plane and the upper (+) half-plane; sc stands for "semi-circle." 
Using this result, v^(Q becomes 



(r) e -rC ( f M^{t-r)e^dt f Ai£> (t; r)e Tt dt 



2 ^ \JCe t-C JC 2 t-C J 

I , e ^(Aig)(t;r)-Ai^(t;r)) , (m £> (t; r) - Aig (t; r)) 



.. . . dt\ (4-13) 

2»tt \Jc 6 t-C Jc 2 t-( 



I , e rt (4'(t;r)-Aig;'(t;r)) ^ | , 
_ e"^ / f Ai£J (t;r)e Tt dt f Ai ( £ (t;T)e Tt dt f Ai { £ (t;r)e Tt dt\ 

It can be noticed that, in all three terms in the last expression, one can close the contours C\,C2, and C@ by 
adding the corresponding arcs of infinite radius. For instance, C\ can be made a closed contour by adding 
the arc {re e \0 € [27r/3, 47r/3]; r — > oo}. Such addition is allowed because the term e^ s+r ' t in 



Aff(t;T) e * _ 1 

t-C 27Ti(t - () 



C.i 



( S + r) r e*( s+r )- s3 ds (4-14) 



is suppressed over the infinite arc when s is integrated over C4. The other terms work the same. Then 
each integration on the closed contour becomes a residue calculation that leads to the definitions in the 
theorem. □ 

4.2 Asymptotic behavior of .4^(0 for large £ and growing r 

The matrix A r (() will be used to build an appropriate local parametrix near z — j3. The resulting global 
parametrix will have a jump across dD c . To control this jump it is necessary to see how A r {() behaves as 
C -> 00 (recall from (3-9) that ( = ©(n 1 / 6 ) on <9D C ). Special care is needed because r is also growing with 
nasr = (9(n 7 ) for < 7 < 1/12. 

Let A( r .ij) stand for the (i,j) entry of A r (C)- We will show in Section 6 that the leading-order asymptotics 
of the kernel are smooth near z — (3. Therefore it is sufficient to restrict ourselves to £ in region /. The 
proofs of the following two propositions are given in Appendix A. 

Proposition 4.2. As ( — > 00 with ( S / and for r = C(n 7 ), < 7 < 1/12, the entries of the first and 



1G 



second columns of A r (() behave as follows: 



6M+5 



/ r 4M+4\ 



A r MC) =d C A${C,T) = ^(C 1 / 2 -r)-C 1/4 e"^ 



.2,3/2 



6M + 5 , , / JM+4\ 

-i+E^ ,a V J/a + o(^r) 



. (4-15) 



(4-16) 



^ i81 (C) = Ag- 1 )(C;r) = ^L(C 1 / a 



2 ,3/2 



AM2)(C) = A#(C;r) = - ^ (C 1/2 +r) r C 1/4 ei 



-3/2 



1 2 6M+5 

J- + — + — + V j4 (r.31) f -3/4-i/2 
£3/4 ^ £5/4 T £7/4 T J ■» 

/- r 4M+4 ^ " 



(^£3M+15/4 J 



£3M+15/4 
6M+5 

1+ ^ A^ 12) (-i/ 2 + 



3 =3 
6M+5 



„4Af+4 



£3M+2 



Diw -f-o /• 4M4-4 \ 

^, 22 )(0=5c^ ) (C;r ) = -^(c 1 / 2 + .) r C 1/4 e ^ 3/2 i + E ^^(j^) 

j— 3 

A. 32 )(C) - 4r X) (C;r) = --^ (C 1 / 2 +r) r e^ 



(4-17) 
, (4-18) 

(4-19) 



1 T T 

£3/4 _ ^5/4 + £774 



2 6M+5 



+ e 4 r,32) c 



(r,32)._3/4-j/2 



J =3 



+0 



r 4M+4 \ 
£3M+15/4 J 



(4-20) 



where the A'f' 1 ^ are polynomial functions of r and t. 



We will not need the exact form of the coefficients Aj , but they can be determined algorithmically 
through the proof in Appendix A. 

Proposition 4.3. The entries of the third column of A r {() behave as follows as ( — > 00 and for r = 0(n 7 ) 7 
< 7 < 1/12; 



(-l) r+1 H c-< 
2i^ 



M 



K1 



M+l 



where 



Vj (r):= £ I AigfrrJe'Vdt. 

f 4 R1 "^ C B 



(4-21) 



(4-22) 



(r) 

The asymptotic expansion of -4( r> 23) = cfe 3 is (/wen in the same form but with different coefficients: 

Vj(r)-^fij{T)-=-TVj(T)-(r + l+j)Vj-i{T), »?_i(t):=0. (4-23) 
TTie asymptotic expansion of A T ^z is given by At r -i,iS) using (4-21). 
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As a side remark, by applying the rotational symmetry of the Airy function, we get (2-8) repeated here: 



%(t) 



i(t) (e 



■t/u 



dt. 



exp 



2m 

IT 



(4-24) 



Evaluating the integral, 770(7-) is 1 at r = 0. 

The first few terms in the leading-order expansion of the entries of A^(() are given in Section A. 3. 

For r = 0(n 7 ), < 7 < 1/12 and £ £ I, Propositions 4.2 and 4.3 yield immediately the following expres- 
sion for A r (C)' 



MO 




l (0) (C)S(C)C(C), 



(4-25) 



where S(() is given by (3-11), 



C(C) := 



V2 e -IC 3/2 











(4-26) 



2iti 



and 



L ( o)(C) = I + 



O 



o 
o 



o 



(4-27) 



There are several important points to note concerning L( )(C), which we will use to define the global 
parametrix inside D c . We are specifically interested in the behavior of L( )(C) f° r C £ <9B C , where £ = ©(n 1 / 6 ). 
The observant reader will note that several of the terms in L( )(C) are not decaying as n — > 00 for ( £ <9B C . 
Ideally we would have L/ )(£) close to the identity for large n. This will be achieved by a series of pseudo- 
Schlesinger transforms (from now on, we simply use the term "Schlcsinger transform"). Each transform 
only impacts one column at a time and will modify L( ) into a new Lyj where the most dominant term in 
that column of the expansion (4-27) is sequentially removed. At each step the rate of convergence for that 
column will be improved. 

We will show in detail how such a transform is used to define a new matrix L(i)(£) with the same entry- 
wise growth as L( )(C) except that the highest-order term in the third column (the explicit — rn 1 / 6 /^ term 
in the (1,3) entry) is removed. 

A note is in order concerning the ^-independent term in the (2,1) entry. There is a nilpotent matrix No 
such that L( )N is the same as L( ) except this term is removed. We then redefine H( ) ■= H( )N and 
L( ) := L( )N as our starting point. 

We then apply a finite sequence of Schlesinger transforms to the improved L^j (() matrix to yield a final 
matrix L( 2 )(C) of the form l + 0(n~ 2 / 3 ). In fact it is possible to use such procedure to change L( )(C) into a 
matrix arbitrarily close to the identity. However, other factors in the calculation (see (6-19)) will introduce 
an 0(n~ 2 / 3 ) error into the kernel computation, so there is no point in removing smaller terms. 

We start with two observations. First off, L(q) (C) has no jump discontinuities, and the expansion of 
I J (o)(C) has only negative integer powers of £. This is important since the Schlesinger transform will remove 
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pole terms. Secondly, L( )(C) has an asymptotic expansion for large n. More exactly, once the pole of the 
lowest order is removed from a given entry, the remaining terms are asymptotically smaller as n — > oo than 
what was removed. Furthermore, each entry can be made to decay at an arbitrarily fast rate by removing 
a finite number of terms from the Laurent expansion of that entry. This follows from the behavior of the 
entries of A r {C,)] see Propositions 4.2 and 4.3. The leading-order entries of L( )(£) can be read off from the 
formulas in Section A. 3, and in principle any term can be computed as described in the proofs. 

Now is an appropriate time to explain where the limitation 7 < 1/12 arises (recall r — C(n 7 ) for 
< 7 < 1/12). Note in Proposition 4.2 the terms (C 1 ^ 2 ± i~) r . To explicitly compute the entries of L( )(C), 
it is necessary to use the expansions 

1± T Y-1 + rT 1 r(r-l) T 2 r(r-l)(r-2)r3 /r 4 \ 

1± -£w) -^^72+ — Y( — ± slc^ VCV' ( } 

which are asymptotic expansions for £ £ <9B C only for 7 < 1/12. Note that the number of terms in L( )(C) 
that need to be removed to leave an error of 0(n -2 / 3 ) will depend on how fast r is growing. Indeed, if 7 is 
close to 1/12, then a very large number of terms will need to be removed! 

5 The error analysis 

5.1 Global Parametrix Construction 

The asymptotic expansion of the outer parametrix *&(z) on the boundary of D c (Lemma 3.3) will guide our 
first attempt at defining the global parametrix . Around all of the band endpoints except for (3 we will use 
standard Airy parametrices which match with *(z) up to 0(l/n). For more details on the Airy parametrix 
in a 3 x 3 Riemann-Hilbert problem sec, for example, Bleher and Kuijlaars [16]. For convenience, we will 
denote the Airy parametrices by P A j(z). We also define B A j to be the union of 2g+ 1 small, fixed-size disks 
centered around a\, f3\, a g , f3 g , and a in which the Airy parametrices will be used. We set 

r*(*) = H (0 )(*)S(0 z^(D c UD Ai ), 
*^)(*):= <^H (0) (z)D(z)- 1 L (0) (C)S(C)D(z) z e D c , (5-1) 
|p Ai (z) z G D Ai , 

where H( )(z) is defined by Lemma 3.3, L( )(C) is defined by (4-25), S(£) is defined by (3-11), and 

no \ 

D(z) := 1 . (5-2) 

\0 e nh ( z )- nh (fi)- T ( I 
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Direct calculation shows that ^f^(z) satisfies exactly the same jumps as W(z) inside '. 



(o)+ 



(o)+ 



(0) + 



(0) + ^ 




(z) 



(*) = * 



(0)- 



4,3/2 

-ez 1 ' 1 



_ e U 3/2 +nh(z) 












(-l) r 

_{z) ( (-ir +1 




1 

I C 3/2 








zer 2 n D c , 



z g r 3 n ] 



z e r 4 n i 



(5-3) 



We now define the error matrix E( )(z) by 



E 



(0) 



(5-4) 



Note that E( ) (z) has no jumps on the contours inside D c . We take the boundary of D c to be oriented 
clockwise (i.e. inside/outside is — /+ respectively). Then, using the expansion (3-12), the jump for E( j(z) 
is given on D c by 

V$(z) := E (0) (z)l 1 E (0) (z) + = *^(z)_*-( z); i = H (0) (z)D(z)- 1 L (0) (C)D(z)H (0) (z)- 1 . (5-5) 

We would like this jump matrix to be close to the identity as n — » oo for ( e D c . However, as noted 
previously, L( )(C) does not decay to the identity for ( = ©(ra 1 / 6 ) (see (4-27)). To remedy this situation, we 
will modify L( )(C) in a series of Schlesinger transforms: 



L (0) (C)^L (1) (C)^L (2) «). 



(5-6) 



The first step will be carried out explicitly to explain the procedure and will remove one representative term. 
The second step will remove the remaining error terms up to 0(n~ 2 / 3 ). There is no need to carry out these 
transforms explicitly as we will show they only affect the subdominant terms in the kernel. We will then use 
L(2)(C) to define a refined global parametrix 'S'^(z). The new error matrix E( 2 )(z) defined using W%Jz) 
will be shown to be close to I. 

Our first transform removes the term -rn 1 / 6 /( in the (13) entry of L( ): define a new matrix 



Then it follows that 



10 — 
L (1) (C):=L (0) (C) |0 1 < 



v0 



(5-7) 



L(i)(C) 



(041) 
(0,21) 



(0,12) 



O 



(0,22) 



L 



(0,23) 



o 



yL(0,31) £(0,32) 



L 



(0,33) 



O 



(5-8) 
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Now there is a unique 3x3 matrix F^j, independent of z, such that 

F \ A ^ /6 ' 

F (l) \ T T , X^, X-1 / . . i 



H (1) (z):=(l+-^)H (0) (z)D(z)- 1 [o 1 | D(z) 



z-[3 



! 

jTih—rC, 

H( )(«) ( 1 



F (l) in. .c.i I a -. „ C 



(5-9) 



vO 1 



is holomorphic for z G D c . The matrix F^) can be computed explicitly, but since its exact form will not be 
significant we use the following proposition to assert its existence. 

Proposition 5.1. Fix a function e(n) so lim^^oo e(n) = 0. Assume that we are given a constant q x q 
two-nilpotent matrix N (i.e. N 2 = 0) , a series of numbers {dj(n)\ — oo < j < k} so that dk ^ and 
dj = 0(e), and a q x q matrix H(z;n) that is locally holomorphic at z = and dctH(z;ra) = 1. We also 
assume that H(z; n) — lim^oo H(z; n) = 0(e) uniformly on a fixed, finite disk around z = 0. Then one can 
uniquely determine k constant (in z) matrices Fi,...,Ffc by requiring that H(z;n) defined below is locally 
holomorphic at z = 0: 

U(z;n):=(l+^ + ... + ^lI(z;n)(l + N £ jj . (5-10) 

In addition, 

F J -(n) = 0(e), j = l,...,k (5-11) 

and 

H(z; n) - lim H(z; n) = 0(e) (5-12) 

n— >oo 

uniformly in the same disk around z — 0, and dct H(z; n) = 1. 
Proof. Denote as follows the expansion at z = 0: 

H(z; n) = H (n) + Hx(n)z + H 2 (n)z 2 + • • • (5-13) 
Let us collect all the terms of the order z~ m from (5-10). 

k-m k /k-(m-j) \ fe_ TO 

Terms of order z- m : ]T F ;+m H ; + ]T ]T F^H, N^- + ^ H,N,i i+m . (5-14) 

i=0 J=l-fc \ (=0 / 1=0 

This must be zero for m = l,...,2fc for H(z;n) to be holomorphic at the origin. Since we have only k 
unknown matrices: F l7 F k , the number of equations must be reduced to k equations. We will show that 
the equations for m = k + 1, 2k are contained in the equations for m = 1, k. 

For m > k the first and the last (summation) terms are absent and we have only the middle term (with double 
summations). The middle term is a linear combination of j ^X^f=o™ Fi+fiHjj N|ra = l,...,fcj (which gives 
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an invertible linear system of equations because d k ^ 0) and, therefore, the set of equations are equivalent 
to 







^F ;+m H ; N form =!,...,&. 



(5-15) 



1=0 



These are also obtained by right multiplication of (5-14) for m = 1, ...,k by N because the second and the 
last terms vanish given that NN — 0. Therefore the vanishing of (5-14) for m = l,...,fc is a sufficient 
condition to solve for Fj's. To solve these k equations consider a big (block) matrix of size qk x qk (where 
q is the size of the matrices that appear in (5-14); q = 3 in our case) obtained by adjoining the Fj's side by 
side as [F l7 F 2 , F k ]. Then we can write the k equations into a single matrix equation as follows. 



[Fi,F 2 ,...,F fc ] 



Ho 
Hi 





Ho 



\ 



•0(e) 



/ 



Hfc_i H fe _ 2 • • • H 

fe-l fe-2 

HiNd,+i, H;Nrf ;+2 , HoNd fe -i + HiNd k , H Nd k 



(5-16) 



1=0 



where the middle term of (5-14) is hidden in 0(e). Because dct H(z; n) = det Ho = 1, the big matrix of size 
qkxqk multiplied on the right of [F 1; F k ] is invertible and, therefore, the solution can be uniquely obtained. 
Also it follows immediately that [Fi, FJ = 0(e). From this, it also follows that H(z; n) — linin^oo H(z; 0) 
is uniformly bounded by 0(e). 

Lastly, to show that dctH(z;n) = 1, we take the determinant of (5-10). 



det H(z; n) = det 
= dct 



Fi 

z 



I+ Fl +...+ 



Ffc 

z k 

F fc 



detH(z;n)det (l|N ^ 



it 



(5-17) 
(5-18) 



The left hand side is holomorphic (at the origin) and therefore it must be 1 to match the right hand side. □ 

We now show how to apply Proposition 5.1 to guarantee the existence of the matrix F^). First, using 
the holomorphicity of Q(z) in D c , we see there are unique functions <i(i J )(n) such that 



nh-T( 



d 



Furthermore, using (2-4) and (2-12), 
We shift z by and choose 

T 

H(z;n) := H (0) (z;n), e(n) : 



1/2 ' 



N := 




(5-19) 



(5-20) 



(5-21) 
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Recall from Lemma 3.3 that the determinant of H( ) is one and that H( )(z; n) — lim„^co H( )(z; n) = O(k). 
Therefore, Proposition 5.1 shows that the matrices F( l m ) exist and 

H (1) (z;n) - lim Hm(z;n) = O (^—) . (5-22) 



l (lHZ)"J- 11111 =u — ^ 

Furthermore, any series expansion of a unimodular matrix can be decomposed into products of the form 
[I + N], where N is nilpotent as we show below in Lemma 5.2 These products can then be dealt with as 
explained in Proposition 5.1. 

Lemma 5.2. Consider a unimodular matrix M(z) with expansion 

M(z) =I + Mi/z + M 2 /z 2 + M 3 /z 3 + ..., M 3 eMat 3x3 (C), TrMi = 0. (5-23) 

Then, for arbitrary K > we can find a finite number of nilpotent matrices Nj and integers kj > 1 such 
that 

M(z) = H (l + ^ Vl + Oiz-*- 1 )) . (5-24) 
Proof. The matrix Mi can be expressed as a linear combination of the following nilpotent matrices: 



(5-25) 



010\ /001\ /ooo\ /ooo 
000 , 000 , 001, 000 

000/ \ooo/ \o / \100 

000\ /000\ /l !0\ /0 
000 , 100 , -1 -1 , 1 1 

010/ \000/ \ 0/ \ -1 -1 

We label these basis elements as Nj with j = 1, . . . , 8 so that 
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for some constants {cj}. Now we may decompose M(z) as 



c,N, 



m w= n i+ — m (*)- ( 5 - 27 ) 



The factor M(z) does not have a z 1 term in its expansion: 

M(z) = I + ^ + • • • (5-28) 
The matrix M.(z) is still unimodular and hence M 2 is traceless and can be decomposed similarly as before 

%)= (n(l+^r)j M(*). (5-29) 
One can iterate this procedure to obtain the decomposition to any arbitrary order. □ 
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We can thus apply Proposition 5.1 to each entry in L(i)(£) in (5-8) to remove all error terms up to 
C(l/n 2 / 3 ). Let T(() be the appropriate transform from L(i)(£) to L( 2 )(C) : 

L (2) (C):=L (1) (C)T(C), CeD c . (5-30) 

The matrix T(£) is chosen so that for z on the boundary of D c (where ( — ©(n 1 / 6 )), 

L (2 )(CW) = 1 + 0(^73) , zedn c . (5-31) 

Let p be the order of the highest-order pole in ( which will need to be removed. Then Proposition 5.1 shows 
there are unique 3x3 matrices F (2,1)1 F( 2 . p ), independent of z, such that 

H( 2) (z) := (l+it (T^r) Hj^zJDfzj^TfODfz) (5-32) 



is holomorphic for z G D c . In addition, 



H 



(2) 



(z;n)- hm H(2) (z; n) = (-^) (5-33) 



and 



F (a , m) - O (^3) , (5-34) 
where this (^(r/n 1 / 3 ) error comes from removing terms in the (12) entry of Lm(Q. Now we can define 

'R(z)H (0) (z)S(C) z^(D c UD Ai ), 
:= I H (2) (z)D(z)- 1 L (0) (C)S(C)D(z) z G D C) (5-35) 

R(«)PAi(«) ^»Ai 

where 

The new parametrix , 3>^(z) has the same jumps as W(z) inside D c . Define the new error matrix by 

E (2) (z):=W(z)(*^(z)) _1 (5-37) 



and note that its jump for z £ <9B C is 

(2) 0) = ^C2)\Z)U{Z) n {a) ^)uyz)n (0) 



V^(z) = H (2) (z)D(z)- 1 L (0) (C)D(z)H (0) (z)- 1 R(z)- 1 



(5-38) 

= H (2) (z)D(z)- 1 L (2) (C)D(z)H (2) (z)- 1 . 
5.2 Global Error Computation 

We have defined the global parametrix after two Schlesinger transforms as (z) in (5-35) and the resulting 
error matrix E( 2 )(z) by (5-37). This matrix satisfies a Riemann-Hilbert problem with jumps across the 
contours pictured in Figure 2 and the boundaries of O^j and D c . 
Explicitly the jumps on these contours are: 
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For z on the contours outside of the disks D^j and D c , the jump of E( 2 ) is 



(5-39) 



Here we define the notation 



V< A > := (Jump of ^J-^Jump of W) 

'10 \ / 1 0\ 

1 - e nV 3 (z) ) e -nVi(z) j g 

,0 1 / \ 1, 



L = <9f2 ut n dfl L , 

'l e nVi(z) e nV 2 (z)\ 

1 

v 1 



1 n an- 1 



aliens n an L , 



/q e n(Rc'Pi(z)+iCT) g\ 
1 





dn+ndni, 



(5-40) 



For z on the boundary of D c , the jump of E( 2 ) is given in formula (5-38). 
For z on the boundary of E>Ai> the j um P of E( 2 ) is 

V$(z) = R(z)H (z)S(C) (P A i^)) _1 RC*)" 1 - 



(5-41) 



• On the bands outside of the disks D^j and D c , V| 2 ^ (z) = ^- Also, 011 contours inside O^i and 
vg>(z) = I. 

The next result follows from the definition of g(z; k) in (1-7). 



Lemma 5.3. For sufficiently small n we have the estimate below, uniformly in z over compact sets bounded 
away from the turning points aj, f3j, a, and (3: 

Q (z; K )=g(z)+0( K ). (5-42) 

Lemma 5.4. In the critical regime, the inner and outer lenses have been chosen so that 

(a) On L outside of the disks D^j and B(/3, S) : The real part of P 3 is negative and bounded away from zero. 

(b) On dCtiens n <9£1l outside of the disks D^j and B(/3, 6): The real part of Pi is positive and bounded away 
from zero. 

(c) On the real axis, away from the bands and outside of the disks O^j and D(/3, 8): The real part of Pi is 
negative and bounded away from zero. 

(d) On the real axis, to the right of f3, and outside of the disk B(/3,<5): The real part of P 2 is negative and 
bounded away from zero. 
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Proof. Statements (b) and (c) follow from the analysis of the Riemann-Hilbert problem for the standard 
orthogonal polynomials (see, for instance, [23]). Statement (d) follows from the definition of the critical 
regime, that PziP) is a unique global maximum for P^ix) in R \ [a,/3] (see Definition 1.2). To prove (a), 
recall that Pz{z) — az — g(z) — l\ + £2, and that Pz{fi) — 0. Along the contour L it is clear that Re [az] is 
decreasing. Likewise 

Re[g(z)}= / log \z - s\p min (s)ds (5-43) 

increases along L, as \z — s\ is increasing along this contour for each s £ supp(p m ; n ). See Lemma 3.1. □ 

We will use the following data about the functions V\ (z), V^iz), an d V 3 {z) to control the jumps of the 
error matrices on the contours outside of the disks. 

Lemma 5.5. For k sufficiently small: 

(a) There exists b > such that on L outside of the disk D c : 

Re[V 3 }< -6n~ 3/4 . (5-44) 

(b) There exists b > such that on dQi ens fl dVL^ outside of the disks D c and D^j.- 

Re[V x ] > 6n" 3 / 4 . (5-45) 

(c) On the real axis, off the bands and outside B(/3,<5) andV)^: The real part ofV\ is negative and bounded 
away from zero. 

(d) There exists b > such that on the real axis, to the right of (3, and outside of the disk D c : 

Re[Vi] < -bn- 3 ^. (5-46) 

(e) There exists a b > such that on the real axis, to the right of (3, and outside of the disk D c : 

Re[P 2 ] < -bn- 3 ^. (5-47) 

Proof. Part (c) follows directly from Lemmas 5.3 and 5.4 (c). To prove (a) we divide the contour into two 
parts: one is the section of L outside of D(/3, 6) and the other is the section of L inside of D(/3, S) \ D c . 
We have choosen L such that Re (V^z) — 3 ^ log(z — /?)) is decreasing along L (sec Lemma 3.1). The proof 
of Lemma 3.1 can be modified to show that Re(Pa(z)) is also decreasing along L. When z ^ D(/3, <5), by 
Lemma 5.3 Re {P${z)— 4p log(z — (3)) converges to Pz{z) for k — > 0, hence by Lemma 5.4 (a) there is a n 
small enough so that Re {V${z)) is decreasing along L and is therefore negative on L\D C , and bounded away 
from zero on L \ D(/?, 5). 

From (2-5) one can see that 

V 3 (z) = -\vi(z)+h(z) + ^ logC = h{f3) + 2 -cl ,2 {z-p) 3 / 2 + -^c l (z-rj) + ^ logC+O {{z /3) 2 ) . (5-48) 
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We start by noting that h(/3)<0 for n sufficiently large so this term will only improve our bound. Next for 
z € L n (B(/3, S) \ O c ), we will show that, for n large enough, the |c^ 2 (z — 0) 3 / 2 term dominates the other 
two. We have 



-TJzCi{z - 0) 



-|ci /2 (z-/3)3/ 2 



3r 



< 



3r f 



2c l/2 n l/3|^_^|l/2 - 2c l/2 n l/12' 



and 



f logC 



9k 



I log CI 



< 



3c 



4c 3/2 I* -013/2 - g c 3 



3/2' 



log(n) +C (^ 1/2 ) 



(5-49) 



(5-50) 



\3/2 



<&n 3 / 4 provided 



-|c?^(^-/3)3/2 

for some constant c > 0. 

To conclude the proof of (a) we note that there is a b > such that Re |c^ 2 (z — ft) 3 
that the segments of L lie in the sector tt/3 < 9 < 5ir/3. 

To prove (b) we divide the contour into two parts: one is the sections of the inner lenses outside of D^j 
and O(0, 5), the other is the sections of the inner lenses inside of B(/3, 5) \ D c . That Re [Pi] is positive and 
bounded away from zero follows from Lemmas 5.3 and 5.4(b). Inside H(ft,5) \ D c we have 

4 



■Pi(z) 



-/>\z-P)^ + 0({z-pf) 



(5-51) 



and the result follows as above. 

To prove (d) and (e) we divide the contour into two parts: one is the interval [ft + <5, oo) the other is 
[ft + n~ 1//2 , ft + 8). That Re [Pi] and Re [P2] are negative and bounded away from zero on the first interval 
follows from Lemmas 5.3 and 5.4(c)-(d). Within the interval [ft + n -1 / 2 , ft + 8), (5-51) gives the result for 
(d), and for (e) we have 

V 2 (z) = h(ft) 2 -c\'\z ft) 3 / 2 + ^a(z -ft) + ^\ og ( + 0(( Z - ft) 2 ) . (5-52) 

Again, h(ft) < for n sufficiently large so this term improves the bound. For z e [ft + n~ x / 2 , ft + 8), we will 
show that for n large enough the — |c^ 2 (z — 0) 3 / 2 term dominates the other two. We have 



J/iCxiz- ft) 



-§c?/ 2 (z-/3)3/ 2 



3r 



< 



3r f 



2c l/2 n l/3| z _ 0)1/2 - 2c l/2 n l/12' 



and 



flogC 



-hl /2 (z- 0)3/2 



9k 



log CI 



< 



3c 



log(n) +0 (^ V2 ) 



4c 3/2 \ z _ 0)3/2 - 8c 3/2 

for some constant c > 0. To conclude the proof of (e) we note that there is a b > such that 
Re 



■\c\'\z-ftfl* 



< -bn- z '\ 

We can now give bounds on the jumps V^' (z) of the error problem. 
Lemma 5.6. In the near-critical regime, for large n, 

(a) Off the boundaries o/D c and Dai- There is a constant b > such that 

Vjj } (z) = 1 + (e- fc " 1/4 ) , z i (0D C U 9B Ai ). 



(5-53) 



(5-54) 



□ 



(5-55) 
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(b) On the boundary o/D c : 



1 



Vjj (z)=I + 0[-^), ze dO c . (5-56) 



(c) On the boundary of Dai • 

Vj^(z) =1 + (-), zedD Ai . (5-57) 



(2) 

Proof. Part (a) follows from equation (5-39), Lemma 5.5, and the boundedness of ^f(2)(z). 

Part (b) follows from (5-38) along with (5-31), (5-33), and the uniform boundedness of D(z) inside D c as 
n — > oo. 

For part (c), first recall from the Schlesinger calculations that Fn m ) = ©(r/n 1 / 2 ) and F( 2 , m ) = ©(r/n 1 / 3 ) 
(as follows from (5-11), (5-22), and (5-34)). Recalling the definition of R(z) in (5-36), we see that 

R(z;n)= 1 + (5-58) 

Now from (5-41) we have 

V^(z) = R(z) ^I + oQ^R^)" 1 , (5-59) 

from which (5-57) follows. □ 
Lemma 5.7. In the critical regime, for n large, 

E (2) (z) = 1 + 

uniformly in z. 



7 2 /3 



Proof. Denote := 9D C U <9B^y and T^y the remaining contours on which the error E^(z) has a jump: 
the inner and outer lenses, and the real axis, outside of the regions D c and D a j- From Lemma 5.6(b)-(c), 

vg?(*) =i + o(-^g), rer c . (5-60) 

Then for n sufficiently large there exists a constant c such that 

IIVJ5 -I|U=(r c) + ||vg } -I|Uo. {ro) < cn-V\ (5-61) 
From Lemma 5.6(a), for n sufficiently large there is a constant c such that 

Hvjy - + ||V ( ( 5 - IlU-cr^) < ce" c ". (5-62) 

The lemma then follows by a standard technique that consists of writing the solution to the Riemann-Hilbert 
problem in terms of a Neumann series involving V"( E ) — I ( see, for instance, [26] Section 7.2 or [27] Section 
3.5). □ 
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6 The kernel near the critical region 

Our main goal now is to obtain the asymptotic form of the kernel K n (x(( x ),y(( y )) (see (1-6)) uniformly 
for Cx,(y m compact subsets. The first observation is that K n (x,y) is smooth in x and y because i) the 
first column of Y has no jump and ii) the second and the third rows of Y -1 have no jump. However, it is 
still possible for the leading-order asymptotics of the kernel to exhibit a Stokes-like phenomenon, namely, a 
discontinuous change with respect to parameters. We now show that the leading asymptotics of K n are also 
smooth. 

Combining (3-2), (5-37), (5-35), and (4-25), we see we can write Y for z € D c as 



We note that A-J -1 and Y x diag[e"f v , e§ v , e% ( v ~ 2az )] are the same up to a holomorphic prcfactor. 
Therefore they have the same jump, which is 



This means that the first column of A r J" 1 and the second and third rows of JA^ 1 have no jump. Since the 
leading term in the asymptotic expansion of the kernel is written in terms of those column and rows, our 
asymptotic expression of the kernel is smooth. Therefore, it is enough to consider, say, £ in region I (see 
Figure 3). 

To arrive at our final expression for the kernel we will need to express A r (( y )~ 1 A r ((x) in a simple form 
involving contour integrals. To do so we take advantage of the machinery of the bilinear concomitant. This 
requires writing A r (() as a constant multiple of a Wronskian matrix. Specifically, in region I, 




(6-1) 




(6-2) 




(6-3) 



where 



/ Ai«(C) 
X r (C):= d c A^\0 
\#Ai£>(C) 




(6-4) 



Eq. (6-3) is obtained from the identity 



a c Aig- 1 )( C ) = Aig 



(O-rAig-^C). 



(6-5) 



From (6-1) and (6-3), we can write 




(6-6) 
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Figure 5: The choice of contours for the r-Airy parametrix in region I. 

We will now express X r _i(Cy) ^\ r — i(CaO i n ^ simple form using the bilinear concomitant. 

6.1 Simplifying Xr-i(Cy)~ 1 Xr-i(Cx) using the bilinear concomitant 

This proposition is a specific example of the more general results on the bilinear concomitant given in [8] . 

Proposition 6.1. Recall the contours C\, Ci, and C^, are defined in Figure 4- Let the dual contours C\, C2, 
and be defined as in Figure 5. Then the entries of Xr-i(Cy) Xr-l(Cx) for i 7^ j are given by 

s 3_ t 3 

Proof. In the proof we will use the same notation as [10] Section 3, and the technical details whose proofs 
we leave out are given there. 
Recall 

A^(C) = ^~ f (t + rY-'e^^dt, 3 = 1,2,3. (6-8) 

The main observation is that, for fixed r, each Kv^T (j — 1,2,3) satisfies the same ordinary differential 
equation of third order which we now derive using integration by parts: 

- sb I i (« + T > re< ""' 3 ) dt = 1 1 (' + « - 12 » <( + T » " + 

= ip-f (r + K-a?){d ( +^)(t + T)'- 1 eM'+ T >-''' 3 dt (6-9) 

2«7T JCj 

= (C(«< + r) + (r - dl - rdD) Aig-^CC). 
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Associated to each of these equations is an adjoint equation [29] obtained by replacing > —d( and 
interchanging the order of multiplications and differentiation operators. Its solutions have the form 

Ai| _1) (C) := J s j^ye^-^ds, j = 1,2,3, (6-10) 

where the dual contours C%, C2, and C3 are defined in Figure 5. The bilinear concomitant [29] for this 
differential equation is a bilinear pairing between the solution space of an equation and its adjoint. For the 
case at hand it admits the following double-integral representation [8]: 

B(^- l \ut l) )=^ i J £ J c ((t + S )(t + r) + ^-C) { -^^e^^dtd S . (6-11) 

Although C appears in (6-11) it can be seen by direct differentiation that the expression (6-11) is independent 
of £. We may follow the steps in the proof of Lemma 3.3 in [8] to show that 2 

sfAig-^Aig" 1 ^ = C$C U (6-12) 



where CjjjCj is the intersection number of the contours Cj and Ci. Introduce the Wronskian Xr(C) of solutions 
to the adjoint equations with entries 

(xAOhrn-i-ir-^r^SiO, l<i,m<3, (6-13) 

and denote 

(C -r -1\ 

F(C) := -r -1 . (6-14) 

V-i oj 

Then, by the definition of the bilinear concomitant in (6-11) and the pairing (6-12), 

= Xr-l(C)F(OXr-l(C) = I. (6-15) 



l<i,j<3 

Therefore % r _ 1 F is the inverse of X r -i- Expanding out the entries of Xr-i an d Xr-i m 

(X r -l(C 3 /)- 1 X,-l(C»)) jj = E (Xr-l(Q) lk (F(Cy)) km (Xr-ltt*)) m j (6"16) 

k' 



gives 



(Xr-iCCvJ-'Xr-i^))^ = ^-j £ j c ((t + s)(t + T) + S 2 -( y ) ^^' e^+^-^dtds. (6-17) 

We are interested only in the case i 7^ j and hence the contours of integration in (6-17) have no intersection; 
this allows us to integrate by parts and obtain in the integrand a harmless denominator (t — s). Then we 



2 The paper [10] contains a wrong sign in front of the intersection pairing. In fact formulas (3.38) and (3.39) should have the 
opposite sign in front of the second term in the respective integrands. Moreover, the wrong overall sign was obtained in using 
Lemma 3.3 to derive the final formula since formula (3.37) was to be the difference of (3.39) minus (3.38) but apparently it was 
miscomputed later as (3.38) minus (3.39). 
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have the straightforward chain of equalities, where only integration by parts is used: 

JcJc, V \s + tY [ (s + tY \J 
= [ I (W«) + 7 J \ ft + ^ e^+C-t-Crftd, 

/■ /■ e ^*-C»" / -r o r ,\ (t + r) r .*-*> , , (6_18) 

= / / +t 2 + s 2 )j '—e^^dtds 

JcJc, t- s V + T S + T J (s + tY 

f f (t + rY s 3 -t 3 e Cxt-(ys 

= / rTT e ^( 9t + 5 ^ dtds 

Jc t Jcj ( s + T ) t- s 
= (Cx-Cy) / / T t : 

Jc, JCj (s + tY t-s 

as desired. □ 

6.2 Proof of Theorem 1 

Proof. From (2-4), it follows that 

H (2) {y)-^ {2) (x) = I + G(x - y) = I + O (^f) 

Recalling E (2) (z) = I + C(l/n 2 / 3 ) (see Lemma 5.7), we see 

r 1 OX^/nVe \-! 
r 1 n-V« H (2) (a:)- 1 E (2) (i,)- 1 x 

10 0/ \ I J 

I n 1 / 6 \ / r 1 \ 
xE (2) (x)H (2) (x) n-V6 r 1 

\ 1 / V 1 o / 



x,yeD c . (6-19) 



(6-20) 

= I + o( - ~M 
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Using this and (6-6), the kernel is given by 

\(V{x)-V(y)) 



K n (x,y) 



• 2 ' 



2iti{x — y) 

e -^(V{x)-V(y)) 

2iri(x — y) 
1 



0, e - nV{ - v \e- n{ y^- a ^ Y{y)~ x Y{x) 



v(v) 



e 2 



2iti(x — y) 



[o,i,i]x(C„)- 1 x(Cx) 



c* - c 



n l/3 



x(C*) 








(6-21) 



1 



1 C„3_ + 3 



1 + 



n 



1/3 



(f + r) r es( s * )+Cx*-CyS 



(27ri) 2 (a; - y) Jc 2+ c 3 J Cl ( s + T Y * ~ « 

where, in the last equality, we have used Proposition 6.1. We now use 

Cx ~ Cy 



dtdsll + O 



Cx Cy 

„l/3 



x - y 



= n 2 / 3 Cl ( K ) (l + G(n- 2 / 3 )) 



(6-22) 



which follows from (2-4). 
We now recall that 



C,(«) = C(«; «) - n'C!(«)(« - 0(k))(1 + 0(* - /?(«))) 
= n 2/3 Cl (0)(z - 0(0) - «0) + 0(n- 2 / 3 ) 



(6-23) 



as long as z — /?(0) = 0(n~ 2 / 3 ), i.e. as long as £ z stays finite for n — > oo. 

This implies that we can replace £ x and £ y in the kernel by the leading approximation in (6-23) without 
changing the error in the kernel. The mismatch from this approximation produces an error of smaller order 
than the one already indicated in (6-21) and yields the overall error term C(n -1 / 3 ) advocated in Theorem 
1. 

By dropping the C(n _2//3 ) in (6-23) and using S := ci$Kn 2 ^ 3 in Theorem 1, 



C:=n 2 / 3 Cl (0)(z-(3)-5 



Cz + s 

n 2 / 3 ci(0)' 



and evaluating K(x,y) at 



(6-24) 



(6-25) 



nsci(0) ' " ' nsci(0) 
we obtain the statement in Theorem 1. 

We conclude the proof noting that C\ may be deformed to C and that Ci + C3 may be deformed to C. □ 

6.3 Connection to the kernel in [1] for nonintersecting Brownian motions 

In conclusion we observe that our expression for the kernel near the critical point is the same that was found 
previously in the literature for nonintersecting Brownian motions: Theorem 0.1 of [1] gives the formula 

a r 1 1 { 7 ib ~ Z eWa3+ib3)+ia{w+ix)+ih{w+cy)dhdadw ' 

4tH J J c J c {ia - ry 
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where C is a contour proceeding from ooe 5,ri / 6 to ooe 7 ™/ 6 and such that — it is above C. To see that our 
formula matches (6-26) we first compute the w-integral to find 

ax,g = 7r^ 7 ^ —7 dbda. (6-27) 

{2mYJ c J c (ta~T) r la + ib 

One makes the changes of variables ib ^ t and — ia i->- s. 

i r r (t + r) r e 3( s3 -* 3 )+c^-c 

(2^ J_ iC j iC (s + tY t 



Kn(( x ,Cy) = 7^2 / / 77^7- 7— dtds - ( 6 - 28 ) 



We conclude by observing that iC ~ C (i.e. C rotated counter-clockwise by 90 degrees) and — iC = C (i.e. C 
rotated clockwise by 90 degrees) (with C,C depicted in Figure 1). Note that the dt contour can be deformed 
(in the finite complex plane) as the integrand lacks a pole in i-space. 

A The asymptotic expansion of *4.( r )(C) 
A.l Proof of Proposition 4.2 

Here we prove the existence of general asymptotic expansions for the first and second columns of A( r ) (C) 
with arbitrarily small errors (i.e. (4-15)~(4-20)). First, observe that (4-16) and (4-17) follow from (4-15), 
while (4-19) and (4-20) follow from (4-18). We then perform a steepest-descent analysis on the original 
integral representation of the generalized Airy function 

^c\Q~^r- I (t + T) r e< t - t3 / 3 dt=^r- [ e^sC+^+C'-^ dt. (A-l) 

3 2i7T J c . 2lTT J Cj 

Let us define Z, T, and T to be scaled versions of (, r, and t, respectively: 

Z ~ ^/3> 7 ' := \T/3> T:= ^T/3- ( A " 2 ) 

Using these new variables, the exponent in (A-l) becomes 

r log (t + t) + (t - j = V - logr + rF(T), F(T) = F(T; {Z, T}) := log(T + T) + ZT - (A-3) 

We define the saddle points T s for s = 1,2,3, as the solutions of 8tF(T) = 0. As solutions to a cubic 
equation, the T s 's can be written explicitly. Let us, however, only write their asymptotic expressions at large 
Z: 

-2 



3m— 1 , ^ 

= -^+E^ (1 " j)/2 +^(^ 



r r 2 



Z(3m+l)/2 ' £(3m+2)/2 



T2 = ^+^-J^ + ^+o(-^) (A-4) 



2Z 2Z 3 / 2 2Z 2 \Z 5 / 2 , 

3m- 1 , , r ; 

= VZ+ rfz^-^ + o' ' 
1 T 2 ( 1 



23m/2 ' _£(3m+l)/2 ' ^(3m+2)/2 
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For large r the coefficients have the asymptotic form = O (7"-7 mod ( 3 )) . This can be seen from the cubic 
structure of the equation drF(T) = 0. It should be mentioned that the above bounds are uniform in T 
inside a finite disk around T = 0. 

Given a saddle point T a , one can expand F(T) by 



F(T) = F(T S ) + 1 -^-(T - T s ) 2 + SF(T), SF(T) = 0({T - T s ) 3 ). 



(A-5) 



The standard steepest descent method gives the leading behavior (in large r) by 
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where the overall sign must be determined from the direction of the contour Cj . 
We find the expansions 

Lemma A.l. For s = 1, 2: 
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J, i>3. 



Furthermore, for large r the coefficients satisfy F$ = O (-p mod(3) ) for k = 0,2, and FjfJ = O (rU-k)mod(3)^ 
for k > 3. 

Based on Lemma A.l, we may obtain the higher order expansion using the general formula 
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where the overall sign is related to the contour involved, and the Sk are the polynomials defined by 

(CO \ CO 

^x^j = l + ^rS k ( Xl ,x 2 ,...,x k )z k . (A-12) 

Note that only the even powered terms contribute after the Gaussian integrals, which all come form the 
formula 

x j dx = r ( . (A-13) 



CO 

r i — X rp3 



2 

We require a lemma to control the growth of the Sk(0, 0, rfa/3\,rfa/4\, . . . , r/fc/fc!) for large r and 
Lemma A. 2. We find for an integer M > t/iat 
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Proof. Using Proposition A.l we see that the leading terms in S k are those with the highest powers of fa 
possible, therefore we have 

fM-l f 

a _ „M J 3 H 

Wl-r M(3!) M-1 4! +•■• 

fM-l t 

s ^= rM Jw^^ + - (A_15) 

j-M+1 

*3M+3 ~ r ( 3 I)M+1 + • • • • 

The terms left off are lower order in r and higher order in \ jZ. An application of the fa, fa, and fa entries 
of Proposition A.l gives the result. □ 

The series in (A-ll) does not converge in most cases, and one must use the original quantity (not the 
series expansion) to estimate the error. We can use that e rF ^ is analytic at T = T x and therefore 
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on a finite disk around Xi of order n 1//3 (the radius of convergence is determined by that of F(T)). The 
right-hand side of this expression follows from Lemma A. 2. Note also that only the even terms will contribute 
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to the Gaussian integration. Then the error is given by 
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where in the last step we used Proposition A.l. Note that 
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We will give here a general argument for the existence of the expansion to arbitrary order of v\ . A 

(r) 

similar argument will give the existence of the expansion to arbitrary order of . 
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Each of the terms in (A-19) have a truncated expansion whose error is dominated by O ( r 4M + 4 /£3M+3~j_ 
To see this we begin with 
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and we have used Proposition A.l. One checks that indeed Rj = O ( j-]?//T ) ■ 
Next we examine 
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We control each of the remaining terms in the exponent of (A-22) separately: 

r(6M+6)/jl-l / / I ,-/3 | + i x r(6M+6)/jl\ 

^»'<l<-"= g >(^/. fl g>r'/')* + o((^-) ). (A-23) 

where we have used Proposition A.l and the identity j/3 — (jmod(3))/3 = L.?/3J- This error will be the 
dominant one (for some choices of j). We have 

r b'/3J+i 

(A-24) 

is largest when jmod(3) = 0. Suppose j = 31, then we have the error arising from the expansion (A-23) is 

r (l+i)\(2M+2)/q \_ 

Noting that 

j(l + 1) - 3Z/2 < 2 7 - 3/2, for I > 1, (A-26) 

we conclude that the dominant error contribution to (A-22) from the expansions (A-23) for 3 < j < 6M + 5 
is 

/ „4M+4 \ 

° {^m+s) • (A-27) 

The remaining terms in (A-19) may be analyzed in a similar way, the error term (A-27) remains the dominant 
one, and our conclusior 
error term of the form 



I >H\(2M+2)/l] ) ■ ( A " 25 ) 



(r) 

one, and our conclusion is that v\ may be written as a finite collection of terms of decreasing order plus an 
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/ r 4M+4\ 



^ J • ( A - 28 ) 

This is small (in large n) for 7 < |. A similar analysis may be carried out for t4^(0 to derive an identical 
error bound. The entries of the second row are then computed directly by taking a derivative in £ of these 
expressions; the entries of the third row are also computed directly by substituting for r = r — 1 in the 
formulas for the first row. This completes the proof that the expansions in (4-15)-(4-20) in Proposition 4.2 
are asymptotic. 

A. 2 Proof of Proposition 4.3 

Here we prove the existence of general asymptotic expansions for the third column of A( r )(C) (see (4-21)). 
To obtain the expansion of -4( r; i3) we use the definition (4-9) in terms of Cauchy transforms. 

(r) e"* ( f Ai£(i)e^ f A^(t)e^dt' 
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The above equality is simply the recombination of contours. Then, from Definition 4.1, we can see that the 
numerator is a total derivative, i.e. Ai^(£)c rt = ^Ai^ _1 ^)e rt ^ . Performing integration by parts, we 
get 
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The first three terms are evaluated at the endpoints of the contours, C^.Cii an d C5. Each contour starts from 
the same point and goes to infinity. By Ai^ — Ai^ + Ai^ = and the decay property of each integrand 
at the corresponding infinity, the first three terms vanish. 

We can perform the integration by parts on the remaining three integrals recursively to finally get 
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We intend to obtain the asymptotic expansion of the above in large C- We define A M such that 
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that is, 
Am 



(A,S) JCb 



(~l) r+1 / 1 r + 1 



(r + M)! i M 
r!M! ^+m+i 



dt, (A-33) 



where (A, B) = (6, 5), (-2, 3), (1, 4). 

We divide the above i-integral (on Cb) into two parts: one for \t\ < L and the other for \t\ > L where L 
is set by 

L = n s , < S < 1/12. (A-34) 



There exists n > such that, for n > n , 
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For r = C(n 7 ), < 7 < 1/12, and \C\ — ©(n 1 / 6 ), there exists n a such that the following crude estimate 
holds for all n> n^: 
(-l)H-i 
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This is true since each term in the left hand side is bounded by 0(|£| M ) for n large enough. 
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Using the above two estimates, we conclude that, for n > n (where n satisfies the above two charac- 
terizations), the contribution of \t\ > L to Am is bounded by 

1 e 3 1 1 
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The second inequality holds if we choose no large enough. The third inequality is obtained by e 3* 3/2 dt < 



Jl° e l* 3/2 ^t 1 / 2 ^ for L large enough. The last estimate is for |£| ~ n 1 / 6 . 



IL ^ 2 

Now we consider \t\ < L. 



An error bound of a finite Taylor expansion for a general analytic function / is given by 
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Let us apply this to our case, f(x) = l/x r+1 and y = —t. We get, for \t\ < L, 
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The last factor is bounded by 2 for n large enough because 
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The contribution of \t\ < L to Am is then bounded by 
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where the constant factor depends only on M and r, not on n. 
From (A-37) and (A-41) we get 
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and we obtain the asymptotic expansion 
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This concludes the proof. 
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A. 3 The leading-order terms in the expansion of A.( r )(C) 

We include here the explicit form of the first few terms of the asymptotic expansions of the entries of A( r ) (() 
to give an idea of their form. Note that more terms would be needed to carry out the Schlesinger calculations 
explicitly, but we do not need the exact formulas for our results. We find: 

-WO = ^J (c 1/2 -r)Y 1/4 e-f< 3/2 x 

. r 2 r 5 \ 1 f r 2 T 3rr\ 1 1 r 4 m ( r 2 \ m ( 'r 3 
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Ar,23)(C) = 5 c Ai^(C) - ^-C- r r!e- T S(T) (-£ + O J J , (A-52) 
Ar,33)(C) = Aig- 1 ^) - O - l)!e- T So(r) (l + O (0) , (A-53) 

where t]o(t) is given by (4-22). 
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